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Abstract 



The standard one-parameter scaling theory predicts that all eigenstates 
in two-dimensional random lattices are weakly localized. We show that this 
claim fails in two-dimensional dipolar Frenkel exciton systems. The linear 
energy dispersion at the top of the exciton band, originating from the long- 
range inter-site coupling of dipolar nature, yields the same size-scaling law for 
the level spacing and the effective disorder seen by the exciton. This finally 
results in the delocalization of those eigenstates in the thermodynamic limit. 
Large scale numerical simulations allow us to perform a detailed multifractal 
analysis and to elucidate the nature of the excitonic eigenstates. 
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One of the most attractive problems in condensed matter physics is the locahzation of 
quasi-particles (electrons, phonons, excitons) in disordered matter. The existence or absence 
of the localization-delocalization transition has been found to strongly depend on the system 
dimensionality. The one-parameter scaling theory of localization states that any nonzero 
disorder causes exponential localization of all eigenstates in one-dimensional (ID) and two- 
dimensional (2D) systems, regardless their energies, while in three-dimensional (3D) systems 
only a rather strong disorder causes the state localization. It is to be noted that there exist 
several exceptions to this rule. In this regard, anomalously weak localization is known to 
occur at the band center in ID and 2D ||^J^ systems with off-diagonal disorder. Moreover, 
correlations in disorder may cause delocalization of states in ID systems even in the 
presence of strong disorder. The long-range inter-site coupling also may act as a driving force 
for delocalization in any dimension. We have recently shown the absence of localization 
in a ID Hamiltonian with a special type of long-range intersite interaction, resulting in a 
specific, non-parabolic quasi-particle energy dispersion. Remarkably, the delocalized states 
belong to one tail of the band 0. 

In this work we report further progress along the lines mentioned above and demon- 
strate that extended states may occur in 2D disordered Frenkel systems, where the exciton 
eigenenergies of the ordered system scale linearly with the wavenumber at the top of the 
exciton band. In such a case, the level spacing decreases on increasing the system size in the 
same manner that the strength of effective disorder seen by the exciton [see Eq. (^) below]. 
Therefore, if the disorder is of perturbative magnitude for a given lattice size, it will remain 
perturbative on increasing the size, consequently allowing delocalized eigenstates. 

We consider a Frenkel exciton Hamiltonian on a regular N' = N x N lattice with diagonal 
disorder: 

n = Y.^n\n){n\ + Y: 4m|n)(m|. (1) 

n mn 

Here, |n) is the state vector of the n-th excited molecule with energy and n = {nx,ny), 
UxjUy being integers {—N/2 < nx,ny < N/2, with even). The intersite dipole-dipole 
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interaction is taken in the form Jnm = J/l^^ — ni|^, where J > is the couphng between 
nearest-neighbor (NN) molecules in the lattice (hereafter we assume that transition dipole 
moments of molecules are perpendicular to the plane of the system and that their magnitudes 
are the same). The joint distribution function of a realization of disorder is the product of 
box functions of width A centered around zero. The quantity A/ J is referred to as degree 
of disorder. 

In the excitonic representation, assuming periodic boundary conditions, the Frenkel 
Hamiltonian takes the form 

n = E^k|K)(K| + E(^^)kk'|K)(K'|, (2a) 

K KK' 

where K = {27t /N){kx, ky) runs over the first Brillouin zone, fc^, fcy being integers ranging 
within the interval —N/2 < kx, ky < N/2. Here i?K is the unperturbed exciton eigenenergy 

= J E A e^"""' (2b) 
and (57i)KK' is the inter- mode coupling matrix 



(^^)KK' = ^E^ne^^''-^')". (2C) 



n 



Hereafter we keep long-range terms in ( pE| ) due to their major role. It can be shown 
that near the extreme points of the band K = and K = tt = (vr, vr) the exciton energy 



spectrum takes a linear and a parabolic form, respectively [10|: 



£;k ~ 9.03J-27rJ|K|, |K| < 1, (3a) 
^ -2.65J + 0.4J|K- 7r|^ |K - 7r| < 1. (3b) 

From this it follows that the energy spacing close to the top of the exciton band behaves as 
A^^^ = A/"^^^^, while in the vicinity of the bottom scales as N^"^ = . 

Depending on the degree of disorder and the lattice size, the operator STi may couple the 
extended excitonic states |K) to each other, thus resulting in their localization. Our task 
now is to calculate the typical fluctuation of this matrix in order to gain insight into the 



magnitude of the exciton inter-mode coupling. The corresponding magnitude of interest is 
'^KK' — (I(^'^)kk'P)5 where the angular brackets indicate the average over the distribution 
Y[nP{^n)- After performing the average one gets 

0"KK' ~ O" = --j=. (4j 

Here a is referred to as effective degree of disorder. As we can see, the typical magnitude 
of the inter-mode coupling scales as Af'^^"^ = N"'^. The most remarkable fact is that a 
decreases on increasing A/" in the same manner as the level spacing close to the top of the 
exciton band. This finding has a dramatic effect on the localization properties of the states 
within this region. Indeed, consider for example a finite lattice of size N x N and the two 
first unperturbed states with K = and K' = (27r/A^, 0). The energy difference between 
them, according to Eq. (pa]), is 6E = Att'^J/N. Take now the degree of the disorder to 
be A <^ Air'^J. Notice that this condition is not very restrictive since in actual systems 
the degree of disorder is expected not to exceed the exciton bandwith (~ 11.68J). Under 
this condition, the strength of the effective disorder a = A/N, governing the exciton state 
mixing and thereof localization, is of perturbative magnitude, namely cr ^ 6E. What is 
most important, it will remain perturbative upon increasing the lattice size because both 
magnitudes scale similarly with A^. Hence, these states will not be mixed by disorder and 
will remain extended over the entire lattice independent of its size. It is clear that the 
same conclusion can be drawn for all the states of the linear spectrum range except the 
degenerate ones. They are mixed by any small amount of disorder. However, since different 
sets of degenerate states are not coupled to each other due to the perturbative nature of 
the effective degree of disorder (a ^ SE), the above conclusion is also valid with regard to 
degenerate states. 

Let us now turn to the parabolic range of the energy spectrum, where the level spacing 
decreases as A/"^^ = A^~^ upon increasing the lattice size, i.e., faster than the effective degree 
of disorder a (the same behavior takes place for both edges of the band obtained within the 
NN approximation, namely taking Jnm = when |n — m| > 1). Now, even if one starts 
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with a perturbative magnitude of A at a fixed lattice size (so that a ^ SE), it becomes 
non-perturbative for larger sizes, resulting finally in localization of those eigenstates. 

The results of numerical diagonalization of the Hamiltonian (|1|), by means of the Lanczos 



method |]TT[, unambiguously confirm our qualitative arguments. To examine the character of 
the exciton eigenfunction (localized or extended) we have calculated the inverse participation 
ratio (IPR) of the uppermost exciton state, according to the standard definition IPR = 
Sn l^!^n|^, where the sum runs over lattice sites and it is assumed that the eigenfunction 
"^un of the z/th eigenstate is normalized to unity. On increasing the lateral size N, the IPR 
scales as A^~^ for delocalized states, spreading uniformly over a 2D system. On the contrary, 
localized states exhibit constant values for different A^. 

The IPR as a function of the lateral size is shown in Fig. |I| when the dipole-dipole 
interaction between all molecules is taken into account as well as within the NN approxima- 
tion. The plots comprise the result of 20 averages over disorder realizations and A = J in 
all cases. The slope of the straight line, being equal to —1.91 when dipole-dipole interaction 
is taken into account, is close to the theoretical value —2. This scaling suggests the fairly 
extended nature of the uppermost exciton state. Notice that there is no scaling of the IPR 
within the NN approximation, in perfect agreement with the well-known results stating that 
those exciton states are localized. Figure § shows an increase of the IPR at a threshold value 
13J, suggesting the occurrence of a smooth delocalization-localization transition. 

A comprehensive way to characterize the spatial distribution of eigenfunctions is the 
computation of the singularity spectrum, as explained in, e.g., Ref. |[12|. If we cover the 
system with (M/L)^ boxes of size (in units of the lattice spacing) and define the nor- 
malized gth moments jjik{<i,5) = fJ'li^)/ J2k' f^l'i^) (where 6 = L/M) of the probability 
distribution /ifc(5) = I]„gboxfc l^nP of finding an exciton in the kth box, we may calculate 
the Lipschitz-Holder exponents 

a{q) = lim^/ife(g,5)ln/ifc(l,5)/ln(5 (5) 

k 

which take into account the scaling of the content of each box with the box size, as well as 



the corresponding value of the singularity spectrum 

f{q) = lim^/ifc(g,(5)ln/ifc(g,5)/ln(5. (6) 

k 

The invariance of the singularity spectrum with the system size for a given value of disor- 
der is usually taken |T^ as a proof of the occurrence of the Anderson transition. Nevertheless 



strong fluctuations of the eigenstates of the system near the transition makes it difficult to 
calculate the /(a) curve for the threshold value of A. We have calculated the singularity 
spectrum for values of the degree of disorder below the threshold. Figure |^ shows the broad- 
ening of the /(a) curve with increasing disorder, indicating that the excitonic eigenfunctions 
at the top of the band becomes progressively more and more localized. This result suggests 
again the occurrence of a smooth delocalization-localization transition, as in Fig. 

In summary, we have shown that the statement of the one-parameter scaling theory 
about the weak localization in two dimensions, i.e., that any amount of disorder results 
in localization of all eigenstates, fails near the top of the exciton band where the quasi- 
particle spectrum scales linearly with the wavenumber |K|. The states lying at such energy 
range are delocahzed at moderate strength of disorder and undergo the continuous Anderson 
transition as the disorder degree increases. In our opinion, the failure of the one-parameter 
scaling theory for the conditions considered in the present work is due to the fact that this 
theory deals only with the size scaling of the energy spacing but pays no attention to the 
subsequent renormalization of the disorder (^). As it follows from our treatment, the latter 
effect plays a major role in localization phenomena, violating the one-parameter scaling and 
thus leading to the impossibility to match our results by this theory. 
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FIGURES 

FIG. 1. Lateral size scaling of the inverse participation ratio of the uppermost exciton eigen- 

function at A = J, obtained by averaging over 20 realizations of the disorder. 

FIG. 2. IPR of the uppermost exciton eigenfunction as a function of the degree of disorder for 
a system of size N x N (shown in the plot). Results comprise 10 realizations of disorder. 

FIG. 3. Singularity spectrum /(a) of the uppermost exciton eigenfunction for various degrees 
of disorder for a system of size 96 x 96. Results comprise 5 realizations of disorder. 
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Figure ^ 




